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Abstract 

Some new global results are given about solutions to the boundary value 
problem for the Euler-Lagrange equations for the Ginzburg-Landau model of a 
one-dimensional superconductor. The main advance is a proof that in some pa- 
rameter range there is a branch of asymmetric solutions connecting the branch 
of symmetric solutions to the normal state. Also, simplified proofs are presented 
for some local bifurcation results of Bolley and Helffer. These proofs require 
no detailed asymptotics for the solutions of the linear equations. Finally, an 
error in Odeh's work on this problem is discussed. 

0.1 Introduction 

In 1950 Ginzburg and Landau |1| proposed a model for the electromagnetic prop- 
erties of a film of superconducting material of width 2d subjected to a tangential 
external magnetic field. Under the assumption that all quantities are functions only 
of the transverse coordinate, they proposed that the electromagnetic properties of 
the superconducting material are described by a pair (0, a) which minimizes the free 
energy functional 

G = i / V(0 2 - 2) + ^ + 20V + 2(a' - hf )dx. 
2d J _ d k 

The functional G is now known as the Ginzburg-Landau energy and provides a mea- 
sure of the difference between normal and superconducting states of the material. 
The variable <fi is the "order parameter" which measures the density of superconduct- 
ing electrons, and a is the magnetic field potential. Also, h is the external magnetic 
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field, and k is the dimensionless constant distinguishing different superconductors. 
So-called "type I" superconductors have < k < 4g while k > 4g for type II super- 
conductors. (But this is really only valid for large d; see 0.) 

The existence of minimizers for the functional G is proved in a standard way, and 
such minimizers satisfy the following Euler-Lagrange boundary value problem: 

0" = K 2 0(0 2 + a 2 - 1) (1) 



a" = cj) 2 a (2) 

with boundary conditions 

4>'{±d) = 0, a'(±d) = h. (3) 

It is not hard to show that the solutions of physical interest are such that <p > 
on [—d, d] , and this is the only kind of solution we will consider in this paper. Also, 
h > 0. Our goal is to determine for what values of h, d, and k the problem has 
solutions, and how many solutions there are in various parameter ranges. 

There are two kinds of solutions of interest, so-called symmetric solutions, where 
4>(x) is an even function of x while a(x) is odd, and asymmetric solutions, where these 
conditions are not satisfied. There is a family of trivial solutions, called "normal 
states" of the form 



4>{x) = 0, a(x) = h(x + c), 



which are obviously symmetric when c = and asymmetric otherwise. In an early 
paper [20| Odeh studied when non-trivial solutions may bifurcate off these normal 
solutions. He concluded that symmetric solutions did bifurcate from the branch of 
normal solutions, but as we shall see just before Lemma ||, his argument had a flaw. 
He also considered whether asymmetric solutions could bifurcate off the normal state, 
but reached no definitive conclusion. 

Subsequently, Bolley and Helffer wrote a series of papers on the problem H, [BL 
10| , []TlJ] , fl2|l and other references cited in JTT|. They gave a quite thorough treatment 
of the local bifurcations which can occur from the normal state, with these results 



summarized in |TT|] and [[12]. Among many results, they gave the correct formulation 
of when and how symmetric solutions bifurcate from the normal state, and did not 
make the error made by Odeh, though they appear not to have noticed the discrepancy 
with his assertions. However some of their proofs are complicated, so we will give 
some simplifications. The proofs below are self-contained, and in particular, we note 
that at least for the results considered below, it is not necessary to use detailed 
asymptotics for the parabolic cylinder functions which solve the relevant bifurcation 
equation. 
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A problem of particular physical interest is whether, as the strength of the mag- 
netic field is lowered, asymmetric solutions bifurcate from the normal state before the 
symmetric solutions. This problem is discussed by Boeck and Chapman and by 
Aftalion and Troy [[|. They relate this question to the formation of vortices in the 
medium, a phenomenon that can not be seen in the one-dimensional model. Accord- 
ing to these authors, if d is neither too small nor too large, and if the asymmetric 
solutions bifurcate first, then interference between the two symmetrically placed so- 
lutions at either edge of the slab can produce a row of vortices down the center of 
the slab. The only rigorous result on this problem is by Bolley and Helffer, who show 
that when d is sufficiently large, it is indeed the case that h a > h s . This is one of the 
results for which we give a simpler proof below. A related result in two dimensions 
with radial symmetry appears in ||. 

After formulating our results we received a new paper of Aftalion and Troy 0, 
who did a thorough numerical study of how the bifurcation curves change with k and 
d. Based on these computations they make a variety of conjectures, some of which 
are related to our work. They conjecture in particular that h a > h s for any («, d) 
such that asymmetric solutions exist. (This has not been proved, though it is well 
accepted by physicists @].) Subsequently Aftalion and Chapman have used methods 
of matched asymptotic expansions to study some of the phenomena found by Aftalion 
and Troy [||[. 

The first rigorous study of the global bifurcation diagram for symmetric solutions 
was by Kwong |]19" |. He proved that for any (k, d) there is a unique curve of symmetric 
solutions, which can be given in the form h = /i(0(O)) for < 0(0) < 1. This curve 
is smooth, and h(l) = 0, h(0) = h s . Hastings, Kwong and Troy studied the nature 
of this curve for large d, showing that it has at least one minimum, followed by at 
least one maximum, if k > a/1/2. This implies that for some values of h there will 
be at least three solutions of the boundary value problem in this range of k and d. 
They also showed that for any fixed k G (0, if d is sufficiently large then for some 
range of h there will be at least two solutions. More recently, Aftalion and Troy 
proved that for sufficiently small nd, there is only one symmetric solution, and there 
are no asymmetric solutions ||. (Numerically it appears that asymmetric solutions 
begin when nd reaches approximately .905 

Up to now, very little has been done concerning the global structure of asymmetric 
solutions (in the parameters K,d,h) or of bifurcations away from the normal states. 
Some initial conjectures were made by Aftalion Q. However a numerical study by 
Seydel |^3[ shows that the picture can be quite complicated. He considers only a 
single configuration, namely d = 2.5, k — 1, and presents essentially the graph in 
Figure 1, in which h is plotted against the value of a at the right-hand end of the 
interval [— d, d]. (Seydel uses a(—d) instead of a(d). There are a number of possible 
"bifurcation curves" which one can draw for this problem. For example, we could 
plot h vs 0(0), as was done for symmetric solutions in [IH]. We elect here to follow 
Seydel and plot a(d) vs h. Either kind of curve gives the important information of 
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Figure 1: The horizontal axis is h and the vertical axis is a(d). The solid curve is the 
branch of symmetric solutions while the dotted curve is the branch of asymmetric 
solutions. The end points of these curves, other than (0,0), are bifurcations from a 
normal state. 



how many solutions there are for a given h.) 

Among the features we see here are the existence of up to seven solutions for a 
given h, and bifurcation of asymmetric solutions from the symmetric branch. It must 
be remembered, though, that asymmetric solutions occur in pairs, and modulo a sym- 
metric reflection, Seydel finds up to two asymmetric solutions and three symmetric 
solutions for fixed h. 

There are two obvious questions to ask concerning the Seydel result. How does 
the picture change as d and k vary, and what are the stability and minimization 
properties of these solutions? 

As stated above, the first of these questions was studied thoroughly in M. From 
their bifurcation diagrams, and results in [12] one can infer results about local stabil- 
ity near bifurcation points. Global minimization was studied by Hastings and Troy 
|T7| . In addition to demonstrating the existence of asymmetric solutions for large d, 



they showed that in some parameter range there are asymmetric solutions but no 
non-trivial symmetric solutions. They also showed that the energy of the asymmet- 
ric solutions can be negative, so that a global minimizer of the Ginzburg-Landau 
functional G must be asymmetric. 

We have done some numerical work to consider the robustness of an asymmetric 
global minimizer as we move into the region where both symmetric and asymmetric 
solutions exist. More precisely, choosing the "Seydel" values k — 1, d — 2.5, so that 
asymmetric solutions exist, we started with h = h a , the asymmetric bifurcation point. 
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Figure 2: The conventions for Figure 1 are also used here. The pairs (k, d) are, 
clockwise from the top left: (.2, 1) , (.35, 2) , (1, 1.25.) , (.5, 2) , (.55, 2) , (1, 2.5) . 



As h is lowered, initially there are only asymmetric solutions, and at points along the 
branch of asymmetric solutions we evaluated the Ginzburg-Landau functional G and 
found it to be negative, so that the solutions must be global minimizers. (This is 
in accord with the result of Hastings and Troy.) Lowering h further we reached the 
region where there are both symmetric and asymmetric solutions. At a decreasing 
sequence of h values, we evaluated G at each of the solutions existing for these values 
of h; up to 5 distinct solutions. We found for a considerable distance down the original 
curve of asymmetric solutions that G takes its minimum on this curve. Thus, the 
minimization property of this asymmetric branch appears to be very robust. 

However, this was a relatively crude examination, and by no means a thorough 
study of the (n, d) parameter space. In this paper we consider only the existence of 
solutions of ([[]) — (|3|) , and not the stability properties of these solutions. 

Turning to the structure of the bifurcation diagram as k and d change, we have 
used the program Auto Jl5| to produce many bifurcation diagrams for different pa- 
rameter values. (Aftalion and Troy also used Auto, independently.) Here are some 
samples: 

From these pictures we see that there are at least the following possibilities: 

1. a single- valued curve of symmetric solutions in the (h, a(d)) plane, and no asym- 
metric solutions, 

2. a single- valued curve of symmetric solutions, from which bifurcates a C-shaped 
curve of asymmetric solutions, 
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3. a C-shaped curve of symmetric solutions, no asymmetric solutions, 

4. a C-shaped curve of symmetric solutions and a C-shaped curve of asymmetric 
solutions, 

5. a C-shaped curve of symmetric solutions, from which there bifurcates a W- 
shaped curve of asymmetric solutions, 

6. an S-shaped curve of symmetric solutions from which there bifurcates a W- 
shaped curve of asymmetric solutions. 

We also see that the asymmetric curves can bifurcate from various parts of the 
C- or S-shaped symmetric curves. Some of these features, and others, are discussed 
in more detail in |J. 

0.2 Statement of Results 

First we consider bifurcations from the normal state. Since the normal state has 
= 0, we rescale by letting 

= at/), 

where a = 4>(—d) and tp(—d) = 1. Then, 

f = K 2 (aV + o 2 -l) (4) 

a" = a 2 tp 2 a (5) 

tp(-d) = l,tp'(±d) =0, a'(±d) = h. (6) 

For a = there is a family of solutions (tp Q , h(x + c)) where h = h(c) is chosen so 
that the linear problem 

^ = K 2 {h 2 {x + c) 2 -l)tp , (7) 

M~d) = l, $,(±d) = 0. (8) 
has a unique positive solution tp Q . From standard linear theory []nj we have 

Lemma 1 For each c there is a unique h = h{c) > such that (0) — @ has a positive 
solution. When it exists, this solution is unique. 
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Thus (£|) — (H) is degenerate at a = 0, in the sense that there is a continuum of 
solutions, because c is arbitrary. 



To remove this degeneracy we reformulate the problem. Consider the equations 
@ — © with initial conditions ip(—d) = l,i/)'(—d) = 0,a(—d) = h(c — d), and 
a'(—d) = h. Then consider ip'(d) and a'(d) as functions of (a, c, h) . We wish to solve 
the equations 



ip'(d)=0,a'{d)-h = (9) 



for (c, h) as functions of a. 
Since 



a'(d) — h = I a 2 vjj (x) a(x)dx 
J-d 



we replace Q9) with 



iJ}'(d)=Q, V (x? a{x)dx = (10) 

We get a solution at a = by setting h = h(c) as given in Lemma 1 and looking for 
values of c such that 

r d 

1(c) ■= / (x + c)if>%dx = 0. (11) 
J-d 

Suppose that (]TTD is satisfied for some c = ci. This gives a solution to fllCf ) 
for a = 0, and this solution can be extended to a > provided that at (a, c, h) = 
(0, Ci,/i(ci) the determinant 

det ( I m 

is nonzero. We will see later (equation (PT8|)) that /(c) = implies that ^^f 1 = 0, 
and standard theory implies that ^gjf^ ^ 0. Hence, a unique branch of solutions 
bifurcates from (0,ci,h(ci)) provided that /'( c i) 0- 

For c = 0, a = hx is an odd function, and this implies that ipo is even. Thus it is 
automatic that 1(0) = 0. Further, for a > we can consider instead of (|J) — (||) the 
problem for symmetric solutions. This is (|4]) — (|5]) on [0, d] with ^'(O) = ip'(d) = 0, 



a(0) = 0, a'(d) = h. Kwong's result in |19| shows that this has a unique solution for 



each a 6 (0, 1), defining h as a function of a. If /'(0) ^ 0, then this is also the unique 
solution in a neighborhood of the bifurcation point (0, ci, fo(ci) to © — @ • Hence a 
unique solution, which is symmetric, bifurcates from (0, 0, h(0)) when I'(0) ^ 0. But 



7 



this condition depends on d. When it is not satisfied we expect a more complicated 
picture. 

In |P Of , Odeh claimed that I'(0) is always positive (for any d). He may have 
thought this was so because he thought that ^7 was an even function of x when 
c = 0, but in fact, this function is neither even nor odd. We have the following result, 
also proved by Bolley||, but with a much longer proof: 

Lemma 2 Suppose that for each positive k, and d, and each c, h is chosen as in 
Lemma 1. Then for sufficiently small nd, I'(0) > 0, while for sufficiently large nd, 
I'(0) < 0. 

The (relatively short) proof will be given in section 3. As a consequence we have 

Theorem 3 For any (k, d) there is a bifurcation of symmetric solutions from the 
normal state. For sufficiently large nd and for sufficiently small nd a unique curve of 
symmetric solutions bifurcates from the normal state. In other words, for sufficiently 
small a, and for some 8 > 0, (§]) — @ has a unique solution with \h — h(0)\ < 5, and 
this solution is symmetric. 

Further, suppose that and («2, (fe) are suc ^ that I'(0) < if (n,d) = 

andl'(0) > if(n,d) = (n 2 ,d 2 ), and assume that (n(t) , d(t)) is a real analytic 
curve C joining and (K 2 ,d 2 ), with k(0) = K\, d(0) = d x and «(1) = k 2 and 

d(l) = d 2 . Then there exists a to G (0, 1) and asymmetric solutions (which are nearly 
symmetric) arbitrarily close to (0, hox) with h near ho, k near K,(to) and d near d(to). 
Here h is the eigenvalue found in Lemma [| for (k, d) = (n(t ), d(to)) and c = 0. 



Remark 4 Note that, unlike j^j, we do not need a transversality assumption, and 
with care we could avoid assuming that the curve is analytic. Note also that Lemma 
|^ ensures that suitable points (ACi,di) and (K2,d 2 ) exist. 

The asymmetric solutions obtained in Theorem 3 are, at least initially, nearly 
symmetric, since they start from c = 0. A different sort of bifurcation of asymmetric 



solutions was obtained by Bolley and Helffer fly] , and independently, with a different 
proof, by Hastings and Troy fl^|. In this case we consider a fixed large nd, and vary 
c, looking for other values of c where 1(c) = 0. The symmetry in the problem means 
we only have to consider positive c. It is obvious that for c > d, 1(c) > 0, so using 
Lemma ^| we have: 

Corollary 5 For sufficiently large nd, there is at least one c > where 1(c) = 0. 

In fact, there is only one such c and asymmetric bifurcation occurs at this point. 
Thus, at this positive c where 1(c) = 0, we have h = h a . The uniqueness of this 



positive c was initially shown by Bolley and Helffer flTTJ, but here we will give a 
simpler, self-contained, proof. 
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Theorem 6 For sufficiently large nd there is exactly one c = C\ > such that 
I(ci) = 0. (By symmetry there is also one negative c with this property.) Further, 
h(ci) > h(0). 



Finally, it is not hard to show that bifurcation does not occur for small koI [ [iT|j . 

Now we turn to more global results. The goal now is to show that bifurcation of 
asymmetric solutions can occur from the interior of the symmetric branch, rather than 
just from normal states, and show that the resulting branch of asymmetric solutions 
can be continued in the (h,a(—d)) plane ( d large and fixed) to the asymmetric 
bifurcation point which was found in Theorem |6|. To state this result we must first 
recall a result of Kwong This result is about symmetric solutions, and concerns 



the global bifurcation curve of symmetric solutions, for any fixed d and k. Since we 
are considering only symmetric solutions, we consider ([]]) with the following initial 
conditions: 

0(0) = a, 0'(O) = 0, a(0) = 0, a'(0) = S, (12) 

where a G (0, 1) and 5 > are to be chosen such that > on [0, d] and (f>'(d) = 0. 
By continuing the resulting solution with <fi even and a odd to the entire interval 
[—d, d] we get a symmetric solution to (P , (@) ■ Kwong's result is: 

Lemma 7 For each a G (0, 1] there is a unique 5 > such that the solution of (|ll,|ll"D 

is positive and satisfies 4>'{d) = 0. 

Hence, for each a we obtain a unique h = a'(d) and a unique a(d). Plotting a(d) 
vs h gives the global bifurcation curve for symmetric solutions in the form referred 



to above. An alternative form as used in |18j is to plot h vs a. The main new result 
of this paper is: 

Theorem 8 // the product nd is sufficiently large, then bifurcation of asymmetric 
solutions occurs somewhere along the curve of symmetric solutions. For any fixed k, 
if d is sufficiently large, there is a continuum of asymmetric solutions which connects 
the curve of symmetric solutions to an asymmetric normal state. 



Remark 9 Here we are identifying points in the (a(d),h) bifurcation diagram cor- 
responding to asymmetric pairs of solutions. The asymmetric normal state referred 
to in this theorem must be the one discussed in Corollary and Theorem since 
to within a reflection there is only one asymmetric bifurcation point from the normal 
state. By real analyticity the continuum in this theorem contains a curve which is 
parametrized (in some sense) by h. 
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Remark 10 Also, it is expected that for k > 4g, if d is sufficiently large then the 
curve of symmetric solutions is S-shaped, so that for some values of h there are three 
solutions. In [18\ it was shown that there are at least three solutions, which is con- 
sistent with this conjecture. Theorem Q shows that bifurcation must occur somewhere 
along this curve, but we are not able to prove anything about the location of the bi- 
furcation point on this curve. Similarly, for k < ^ there are at least two solutions 
for large d, and the bifurcation seems to occur on either of the two branches. 



Remark 11 A problem which we have not been able to solve is to determine the 
direction of the bifurcation. As a result, we have not been able to prove that there 
are some values of the parameters (d, k, h) where there are five distinct solutions, as 



was seen in Seydel's original numerical result. In [Tj\] there is a discussion of the 
stability of solutions bifurcating from the normal states. This involves consideration 
of the energy functional G, and we have not studied this topic here. 



1 Proofs 

1.1 Local Results 

In this section, we will denote ipo by if), since we are only considering the linear 
equations. Thus, if) is assumed to satisfy © — © • 

Proof of Lemma ^| 



In ([7]) — (||) we can make the change of variables x — > kx, which removes the term 
k 2 from the differential equation, while replacing d with D = nd, and also introducing 
a new h and c. The latter changes are immaterial in this and subsequent results about 
bifurcation from the normal state, so without loss of generality we will simply assume 
in (|7|) — (|8|) that k — 1, and to remind us of this, replace d with D. Multiply (|7|) by 
if)' and integrate by parts to get 

2h 2 I(c) = (h 2 (D + c) 2 - 1) ip(D) 2 - (h 2 (-D + c) 2 - 1) ip(—D) 2 . (13) 

Consider (|7|) with the initial conditions if>(—D) = l,ip'(—D) = and denote the 
solution by if)(x, h, c). Let P = §f and q = Then 

p" = (h 2 (x + c) 2 -l)p+2h(x + c) 2 ip, p(-D) = p'(-D) = 0, (14) 

and 

q" = (h 2 (x + c) 2 -l)q + 2h 2 (x + c)ip, q(-D) = q'{-D) = 0. (15) 
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Lemma |1] tells us that (|7|) — (|8|) define h as a function of c. We can also see this 
locally by applying the implicit function theorem, solving the equation 



if)'(D,h,c) = 

for h as a function of c. We can do this if p'(D) ^ 0. Multiply by if) and ([7|) by 
p and integrate from — D to .D. With the boundary and initial conditions for if) and 
q, this gives 

if)(D)p'(D) = 2h (x + c) 2 ip(x)dx > 0. (16) 

J-D 



Hence h is a smooth function of c. Further, 

~dc ~ ~p'(D)' 



(17) 



Now multiply (|15|) by ^, (0) by g, subtract and integrate, and use the boundary 
conditions again to get 



r D 

if)(D)q'(D) = 2h 2 / (x + c)ip{x) 2 dx = 0. (li 



This shows that ^ = whenever 7 = 0. In particular, this is true at c = 0. 

Now differentiate (|l|) with respect to c and then set c = 0. Since, then, if>(±D) = 
1, we get: 

h 2 I'(0) = 2h 2 D + (h 2 D 2 -l)q(D). (19) 
From (0) and (fH| ) the equation obtained at c = is 

a r n 2h 2 sif){s) 2 ds 

(ft- «J? ■ (20) 

Since ^ is an even function at c = 0, the right side of ( pOD is zero when x = D. 
Further, the integrand is negative for s < and positive for s > 0, and this implies 
that the integral is strictly negative for —D<x<D. 

Now return to (0) . We see that if)" > —if), and since if)'(0) = 0, this leads 
to if)(x) > if)(0)cosx on [0, .D] for small D. It is therefore seen that as D — ► 0, 
max^gfo^] |V ; ( ;;c ) — 1| — 0. Integrating the right side of the differential equation for if) 

then shows that lmiD^o /iD = so fa — > oo. From (p0|) it follows that = O(h) 

as D — > 0, so g(-D) = O(hD) = 0(1). Hence from fllTf ) we see that for sufficiently 
small .D, J'(0) > 0. 

To complete the proof of Lemma there remains to show that I'(0) < when D 
is sufficiently large. We first need to show that h is bounded as D — > oo. In fact, it 



11 



is well known, e.g. |T7|], that ft — > 1 as D — > oo, but the following technique quickly 



shows that ft is at least bounded: Let p = Then 

p> = h 2 x 2 - 1 - p 2 , p(0) = p(D) = 0, (21) 

with p < in (0, D). (Remember that we are only considering c = here.) It is easy 
to see from fl2~T|) that if /z. — > oo, then D — > 0. Hence, as D ^ oo, ft must remain 
bounded. 

It is also well known that ft > 1, but for completeness here is a quick proof: 



Compare p from (|2l|) with the solution a = — x of 

a' = x 2 - 1 - a 2 , cr(0) = 0. 

An easy comparison shows that if ft < 1 then p(0) = er(0) = implies that p < a for 
all x > 0, which contradicts p{D) = 0. 

Lemma ^] now follows from ( p~9| ) if we can show that q(D) does not tend to zero 
as D tends to infinity. To do this we again use (p0|) , and also (|7|). Multiply (|7|) by ip' 
and integrate from — D to x. This, with fl8|) and ( ^0|) , leads to 

i)'<(*V-l). 

Since ft is bounded for large D, there is some interval around x = of fixed length 
„ > in which ( J)' < -I. In addition, i is decreasing on the entire interval [- A fl], 

so ^ < -|p. But if>(D) = 1 when c = 0. Then (JT^) shows that J'(0) < for large 
D. This proves Lemma 0. 

Proof of Theorem |3| 



Note that I'(0) is a function of k and ci. (We are no longer assuming that k — 1.) 
In fact, it is a real analytic function of k and ci, which is seen by observing that one 
can use the implicit function theorem to prove that ip is a real analytic function of 
(n,d). [14| Thus I'(0) is a real analytic function of t along the curve C, since C is 
real analytic. Since /'(0) does not vanish identically, being nonzero at the endpoints 
of C, its zeros are isolated. Therefore there is a t G (0, 1) such that I'(0) has a strict 
change of sign as we cross (n(to), d(t )) on C. 

We now look for solutions (p = a(ip + w), a = h(x + c) + p , where tp is the 
eigenfunction at t = to, ft is close to ho, and the numbers a and c and function w and 
p are small, and where w is orthogonal to ipo and p is orthogonal to 1, over (—d,d). 
Equation flU) can be written as 

-(V>o + w)" = « 2 (^o + w)(l - a 2 - a 2 (^o + w) 2 )- (22) 
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We now use Fredholm alternative theory and the implicit function theorem in a 
standard way ( ||22|| , chapter 9) to show that for every (a, c, h) with a, c, and {h — h ) 
sufficiently small, there is a unique pair (ft, 7) such that the equations 

-(iPo + w)" = n 2 (ij + w)(l- (h(x + c) + pf - a 2 (ij + wf) + ftifj (23) 
p" = ct 2 U)Q + w )\h(x + c)+p)+'i (24) 

have a unique solution (w,p) such that w'(±d) = p'(±d) = 0. In other words, we 
subtract off suitable multiples ftip of ipo an d 7 • 1 of 1 in equations (|23D — (p4j) 
respectively in order to find solutions (w,p) in the space 

{io G C 1 !— d, d] :w is orthogonal to ip ,w'(±d) = 0} x {p G C 1 !— d, d] (25) 
and p is orthogonal to l,p'(±d) = 0} 

We obtain w, p, (3, 7 as smooth functions of a, c, h. Here we are really solving a pro- 
jected equation. This is a Lyapunov-Schmidt reduction. We want to find solutions 



of (H) - (|D with = 7 = 0. 
Let p denote the solution of 

p" = V> 2 x, p'(±d) = (26) 

with mean value zero. Using (fj) it is easily shown that 

p(a, h, c) = ha 2 (p + o(l)) 

as (a, /i, c) — > (0, /io, 0). Note here that sip 2 (s) has mean value zero. Hence, 

a 2 = (h(x + c) + a 2 (p + o(l))) 2 = h 2 (x + c) 2 + 2a 2 h(x + c)p + terms of higher order. 

(That is, higher order in a.) Note that as (a, h, c) — > (0, ho, 0), w and p tend to zero 
uniformly in [— d, d]. 

Now multiply ( |23| ) by ^ and integrate, using the orthogonality of w and ip , to 

get 

d rd rd 



{h-K){\ s 2 M s )ds + o(l))ds + / o(l)ds = p M s ) ds , ( 27 ) 

J-d J-d J-d 

where the o(l) terms are terms in u> and p which are smooth and tend to zero as 
(a, h, c) — > (0, ho, 0). (Also, the derivative with respect to h of the second o(l) term 
is small if a is small, h is near /to and c is small.) To obtain a solution of ( [22] ) we set 
ft = in (p7D and use the implicit function theorem to solve this equation for ft as 
a function of (a, c) near h = ho, a = c = . The coefficient of (h 2 — h^) in the first 
term on the left of (p7|)is of order 1. For each small a^0 and each small c there will 
be a unique h = h(a,c) near to h for which (|27|) implies that /3 = 0. 
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Now integrate (0) over [—d, d] and use the boundary conditions to obtain, upon 
dividing by a 2 h the equation 



R(a 1 h 1 c) = ^— / 4> 2 a dx = I ((x + c)iPq(x) + o(l))dx = 




-d J -d 



where the small o(l) term is a smooth function of a,h,c and is zero if a = c = 
h — ho = 0. Note that R(a, h, c) = is equivalent to 7 = 0. Also, R(0, ho, c) = 1(c). 
Letting R' 3 (a, h, c) = we have for large nd, by Lemma 0, that .R^O, ho, 0) < 0, so 
R' 3 (a, h, c) < if a and c are small and h is close to h Q . Also, R(a, h, 0) = because 
then both iu and ^0 are even functions of x and a is an odd function of x. Hence, if 
I'(0) < 0, then R(a, h, c) < if a is small, c is positive and small and h is near to ho- 
Similarly, if I'(0) > 0, then R(a, h, c) > if a is small, c is small and positive, and h 
is near to h . 

We now solve fl27|) for ft, as a function of c and a, for small positive a and c. (We 
do not know the sign of h — ho-) We do this for (k, d) on the curve C with (k, rf) 
close to (n(to), d(to))- On this curve, on one side of (n(to), d(to)), I'(0) < and hence 
R(a, h(a, c), c) < if a and c are small and c > 0. How small a and c must be 
depends on the particular point on the curve C. For a given t\ < t but close to t , 
find positive a\ and Ci such that R(a, h(a, c), c) < if < a < ai and < c < c\. 
Choose a fixed t 2 > t but close to t , and then lower a.\ and Ci if necessary so that 
R(a, h(a, c), c) > if < a < a\ and < c < ci.Then somewhere between t\ and t 2 , 
as we keep a and c nonzero and fixed in (0, ai] and (0, ci] and move along C, R must 
equal 0, which gives the required solution. (Note that ho varies continuously with 
(k, d) but this does not affect the argument since everything varies continuously along 
C. ) Since we can choose c arbitrarily small, the solutions will be nearly symmetric. 

Proof of Theorem ||. 

As previously, for the linearized problem (0) — (|TTD we can rescale to eliminate k, 
so we will again assume that k = 1 and replace d with D. We saw in the proof of 
Lemma that for large D, 1(0) = 0, I'(0) < 0, and that consequently Corollary [| 
holds. Further, the definition of 1(c) shows that 1(c) > for c > D. 

From (0) — (H) , with ip > 0, we see that r(x) = h(x + c) 2 — 1 must change sign in 
(—D,D), and from ( |13D it follows that if 1(c) = 0, then r(—D) and r(D) must have 
the same sign, so r() has exactly two zeros in (—D,D), with ip"(x) changing from 
positive to negative and back to positive as x increases from — D to D. Therefore, ip 
has a local maximum at some x G (— D, D), with ip' > on (— D, x ) and ip' < on 



(xoP). 



Lemma 12 For any c G (0, D), h > 1. 
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Proof: Let p(x) = Then p(-D) = p{x ) = p(D) = 0. Also, 

p> = h 2 (x + c) 2 - 1 - p 2 . (28) 
Further, let o~(x) = —x — c. Then 

a> = (x + c) 2 - 1 - a 2 , cr(-c) = 0. (29) 

But p(—D) < a(—D), and an easy comparison of fl28|) and (Fj9|) shows that if 
h < 1, then p < o~ on [— D, D], so p(-D) < 0, a contradiction. This proves Lemma 12. 

Lemma 13 For sufficiently large D, (|7|) — (TTTJ) /ias no solution with ijj > and 
0<c<f. 

Proof: Recall that 7(0) = 0, I'(0) < 0. Let c a = inf{c > 0|/(c) = 0}. Then 
I'( c i) > 0. Assume that c\ < y. We will obtain a contradiction by differentiating 
(|l3r) with respect to c and setting c = C\. 

Recall that 4| = whenever /(c) = 0. Therefore, from (|T3| ) we get 

2h 2 I'( Cl ) =2h 2 (D + cMD) 2 

+2(h 2 (D + ci) 2 - l)ij>{p)q(D) + 2/i 2 (D - c^-D) 2 

where q = |p, so q satisfies ( pT5| ) . From ([15]) and (0) — ([8]) we obtain (f2^) with s 

replaced by (s + Ci), and as in the proof of Lemma fj, we then see that ^\ < 

in (—D,D) and < /i 2 (x + Ci) 2 — 1. Also, as in Lemma fj it is seen that h is 

bounded as D — > oo, and this leads to a negative upper bound of the form 

< -V < 0, (31) 



where r] is independent of D and ci G [0, y] . 
From J(ci) = and (|T3|) we obtain 



^(D + d) 2 - r v y ~ 2h 2 D 2 - i 

io Vinnnr 



Hence, is bounded as D — ► oo. Then fl3T|) and (|30D show that I'( c i) < f° r 

sufficiently large -D. This contradiction proves Lemma O. 



Continuing with the proof of Theorem [|, we now assume that c > y. The result 
will follow if we can show that I'( c i) > f° r an y solution of (0) — (|11]) in this range of 
c. As before, we assume that (ci, /i, t/>) is a solution, and that "0 has its maximum over 
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[— D, D] at x G (-D, D). Also, as before, we know that |x + Ci| < 1, so x < — j + 
We translate the origin to xo, letting if)(x) = x( x ~ x o) — x{y)i so that 

x"(y) = (h 2 (y + x + cif - i)x(y)- 
Now let p(y) = (a shift from the previous p) and let oj(y) = —p(—y). Then 
p> = h\y + x + Cl ) 2 - 1 - p\ p(0) = p(D -x ) = (32) 

and 

J = h 2 (y-x -c 1 ) 2 -l-uj 2 , w{0) =w{D + x ) = 0. (33) 

It is important to recall that Xq < + 1. 

We now need estimates on p and u, which we obtain from the following result: 

Lemma 14 Suppose that for some constants S and A, with h\S\ < 1 and A large, 
r](-) solves the boundary value problem 

r/(y) = h 2 (y + 5) 2 - 1 - r] 2 , 77(0) = 77(A) = 0, 77 < in (0, A). (34) 

Then, 

V (y) >-hy-/3 (35) 



on [0, A], where [3 = \/h — 1 + h 2 6 2 , and 

n{y) < -h(y + 6) + 1 (36) 

on [0, A - 1]. 
Proof: 

Inequality (j3~5P follows by assuming equality at some y\ G (0, A) and using fl34|) 
to show that r/'(?/i) < — /i. This would imply that r/(y) < —hy — (3 for y > y lj so 
that could not vanish at A. Next, observe that fl3~6|) holds as long as h(y + 5) < 1, 
since 77 < on (0, A). If equality holds at some y\ > 0, then h(yi + 5) > 1 and 
V'(yi) = -2 + 2/i(yi + 5) > 0. Also, 77 < and 77' > imply that 77" > 0. Hence, 
77' > — 2 + 2/1(7/1 + 5) > as long as 77 < 0, and if 77 < on [yi, yi + 1] , then 

+ !) > »7(yi) - 2 + 2 %i + 5) = -1 + /i(t/i + 5) > 0. 

This contradiction shows that 77 = somewhere in [7/1,7/1 + 1] . If 7/1 < A — 1, then 
77 = before y = A, a contradiction which proves Lemma |T1[ 

Applying (^) to p with 5 = x + ci and A = D — x shows that 

V>(£>) = x(.D - x ) < ij(x )e-^ D+cl - 1)2+ ^ X0+Cl)2+D - x °- 1 , 
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while applying fl3"5|) to uo with S = — x — C\ gives 

ip(—D) > 4){x )e-^ D+xo)2 - p(D+xo) . 
Combining these and noticing that 2Dxq < —1x\ and (3 < h, we find that for c\ > D 



2 ' 



< e~ rhU (37) 



for some r which is independent of ci, h, and D. 

We now return to (ISO) . Since ip"(—D) > 0, we have h 2 (D - c x ) > h. Further, 
the proof (following equation (|20| ) that when c = 0, h is bounded as D — ► oo easily 
extends to c > 0, since Xq < — c + 1. From ( |3"7D it follows that a bound of the form 

g(£>) > -Lh m D n 4>(D) (38) 

for some L > independent of C\, h, or D will imply that I'( c i) > f° r large D. There 
are two cases to consider, namely, —c\ — 1 < x < —c\ and —c\ < x < —C\ + 1. We 
consider the first, the two cases being similar. Repeating the derivation of ( P0| ) for 
Ci > we obtain that 



■ [D] -2ti 2 (l +1 +1 — ^— I (s + cMsfdsdx). (39) 



D Jx J-c 4 , { x ) J~D 

In the first of the three integrals with respect to x, —D < s < x < Xq, so ip{s) <if)(x), 
and this term contributes less than 0(D 3 ) to as D — > oo. In the third of the 
three integrals, we use the fact that I(ci) = to write 

D -| fx pD j /"D 

s + Ci)ip(s) 2 dsdx — — I / (s + Ci)vp(s) 2 dsdx, (40) 



Cl 



^(x) 2 y_ Cl v(x 



and because — ci > Xo we again have ip(s) < ip(x) and get a contribution 0(D 3 ). The 
second term in ( |39"D is bounded by 



^(x ) 5 



to 



^(x + l) 2 



/x 
(s + Ci)dsdx. 
D 



Again we let p = ^^-xo) anc ^ no ^ e that P' — ~~ 1 ~~ P 2 > P(0) = 0- This gives a bound 
^fio+i) — e<5 ^ or some ^ independent of D, h, or ci, and so the second term in (|39|) is 
0(L> 2 ) as D -> oo. 

This proves the desired bound ([38|) . Hence ([37]) shows that the dominant term in 
fl30|) is the last one, proving that I'(ci) > for large D, if D is sufficiently large and 
I(ci) = 0. Hence there is a unique C\ > with I(ci) = 0. It follows from fll~6| ) — (|l"8"D 
that /i(ci) > h(0), since 7(c) < for < c < c^. This completes the proof of Theorem 
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1.2 Proofs of Global Results 



We now turn to the global results, about bifurcation from the curve of symmetric 
solutions. We consider the full problem ([[]) — (|^) . The symmetric problem can be 
studied on the interval [0, d]. Let (<p(x, a, 5), a(x, a, 5)) denote the solution of flU) — (0) 
which satisfies the initial conditions 0(0) = a, <f)'(0) = 0, a(0) = 0, a'(0) = 5. Recall 
that Kwong proved in |19| that for each a G (0, 1] there is a unique 5 = 5 (a) such 
that is positive on [0, d] and <p'(d) = 0. 
Now let 

s(x) = ^ and z(x) ^ 



85 vy 35 ' 

Let 5 = 5 (a) and let (0, a) be the corresponding solution on [Q,d\. 

Lemma 15 s'(d) > 0. 



Proof of Lemma 15: 



The pair (s, z) satisfies the system 



9" = K 2 [((f) 2 + a 2 - 1)6 + 2a0/i + 2<j) 2 6] 

fj," = 2 /i + 2a<f)6 ( } 

and 

s (0) = s'(0) = 0, ^(0) = 0, z'(0) = 1. 

It is clear that z, z', z" are all positive on (0, d] as long as s > 0, since a > 0, > 0. 
Also, s"(0) = s"'(0) = 0, while s""(0) > 0, so initially, s is positive. Suppose that 
s(xo) = at some first s > in (0, d], with s > on (0, x ). Multiply (|4l|) by and 
QI|) by s, subtract and integrate. We conclude that 

0s' - S0V = / °(2a0 2 ^ + 20 3 s)dx > (42) 
Jo 

and applying the boundary conditions shows that s'(xq) > 0, a contradiction. Hence 
s > on (0, d] and then from (pE2|) we find that s'(d) > 0, as desired. This proves 
Lemma |15[ 

To continue our study of bifurcation from the symmetric branch, now consider ([!]) 
with initial conditions which are possibly asymmetric, namely 

0(0) = a, 0'(O) = 0, a(0) = 7, a'(0) = 5. (43) 

Let 
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q(x) 


d<t> , , 
= W r(x) 


d<f> . , 

= w s(x) 


d<p 


da 


" m 


da 


v(x) 


da , . 

= w w(x) 


da 

= < x) 


da 


da 


= ~d5' 



p(x) 
u(x) 



Then the pairs (p,u), (q,v) , (r,w), and (s,z) all satisfy the system (pEip . The initial 
conditions are 



(p,p>,u,u>) = (1,0,0,0) 
(g, <?>,<) = (0,1, 0,0) 
(r,r', w,w') = (0,0,1,0) 
(s 7 s',z, z') = (0,0,0,l) 



(44) 



at x = 0. (Thus, these four pairs form a fundamental solution for fl4T| ) 
We wish to solve the three equations 



where 



F(a, (3, 7 , 5) = G(a, f3, 7 , 5) = H(a, /?, j, 6) = 



F(a,/3, 7 ,5) = 0'(d) 
H(a,f3,j,6) = a'(d) -a'(-d). 



(45) 



(46) 



For each a e (0, 1], (a, 0, 0, 5o(a) ) is a solution. Starting at a = 1, where 5 = 0, this 
solution continues uniquely as the smooth curve of symmetric solutions as a decreases 
so long as J ^ 0, where 



J = det 



q'(d) r'(d) s'(d) 

q'(-d) r'i-d) s'(-d) 

v'(d)-v'(-d) w'(d) -w'(-d) z'(d)-z'(-d) 



where in ( |4"ID , (</>, a) is the solution at (a, 0, 0, 60(a)). 

Since is even and a is odd, the initial conditions fl44]) imply that p, u, to, and 
s are even functions of x, while u,q,r, and z are odd functions of x. As a result, J 
simplifies to 



J = As'(d)(q'(d)w'(d) - r'(d)v'(d)). 



Lemma 16 For sufficiently large nd, J changes sign between a = and a — 1. 
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Proof: 

Lemma |l^ showed that s'(d) > 0. Therefore, to obtain a bifurcation point on 
the curve of symmetric solutions, we must show that M(a) = q'(d)w'(d) — r'(d)v'(d) 
changes sign along this curve. 

For a = 1 we have <p = I, a = 0. The equations ( f£TD and (|44| ) can then be 
solved explicitly to show that q'(d) > 0, v'(d) = 0, r'(d) = 0, and w'(d) > so that 
M(l) > 0. 

We now analyze M(a) for small a. First, with (q, v) substituted for (6,/j,) and aip 
for in ([IT]) , we multiply the equation for q" by if), (|]) by q, subtract, integrate, and 
use (U) and (|44]) to obtain 



if)(d)q'(d) = 1 + 2/t 2 a / (a^ u + if) aq)dx. 

Jo 

As a — >■ 0, — > ^o, where 

< = k 2 (^ 2 - 1)^ , Vo(0) = 1, Vo(±0 = (47) 

and ho is the unique positive number such that ( pETf ) has a positive solution. (Earlier 
we referred to the solution ip Q of (pE7|) as t/>, but now we must distinguish ^ from 
■0 = - for a 7^ 0.) Since 

v" = a 2 ip 2 v + 2aaif>q, 

( fP| ) implies that t> — > on [0, <i] as a — > 0, so we consider, instead, -, and see that 
as a — > 0, 

v'id) f d 

► / 2h o xif) qodx, 

a Jo 

where 

<tf = K 2 (/i 2 x 2 - l)g , go(0) = 0, q' (0) = 1 (48) 

We proceed in the same way with (r,w), to obtain, finally, that 

Mia) 1 f d f d f d 

lim — = ( / ipQ + 2hoxijjoRodx — AK 2 hl / xipndx ■ / xipoqodx). (49) 

a^o a 2 O [d) Jo Jo Jo 

In this expression the only term we have not defined is Rq, which is lim^o ~ an d 
satisfies 

R' ' = n 2 (h 2 x 2 -l)R + 2n 2 h xip , R (0) = R' (0) = 0. (50) 

To prove that J changes sign it is again convenient to rescale, letting if)o(x) = 
g(nx), so that 

g" = (AV - l)g, g(0) = I, g'(0) = 0, g'(D) = (51) 
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where D = nd and A = It was shown in the proof of Lemma 2 that A > 1. 

Making the same change of variables in (|49|) , we must show that for large D = nd, 



1 (g(y) 2 + 2^yg(y)P(y))dy<A^ yg{yfdy ■ yg(y)Q(y)dy (52) 



where Rq{x) = P(nx), and qo(x) = Q(kx). Hence, we have 

Q" = (AV-1)Q, Q(0) = 0, g'(0) = -, (53) 

K 

and 

P" = (\ 2 y 2 -l)P + 2\yg(y), P(0) = P'(0) = (54) 



as well as (|5"ID . Multiplying ( |53"D by (7 and (|51~D by Q an d subtracting and integrating 
gives 

qX 1 



9 J K9(y¥ 

and similarly from (|54l) we obtain 

p\> 1 ry 



From these we find that 



1) =2X JWJ S9{s) ds - 



P (v)=g(v) !\hJtp^ dXi (55 ) 
Jo Kg(x) 2 

and 

With these substitutions, (|52|) becomes 

g{y) 2 dy<A\ 2 yg(y) 2 tg{t) 2 dtdsdy. (56) 

JO JO ^V S J Js 

It is tempting to approach this result by studying the asymptotic behavior of g(y) 

as D — > 00. In fact, one can show that ^(7/) — > e 2 point-wise, and further effort 
can refine this result. It turns out to be a mistake, however, to study the result of 

y 2 

substituting e~~ for g in fl56]) , because this function does not satisfy the boundary 
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conditions, and this turns out to make the required estimates much more difficult, or 
indeed, impossible. 

Instead, we proceed directly, and this is possible primarily because using fl5T|) we 
can evaluate integrals of the form J yg(y) 2 dy. In particular, multiplying (|51~1) by g' 
and integrating by parts, we find that 



A 2 I yg{y) 2 dy={\ 2 y 2 -l) 



g{y? g\yf 



2 2 

Substituting this in (|)6|) and using the boundary conditions gives 

,9 H 1 f° ,,■>, , X 2 D 2 -1 . N9 f y 1 , lg'(y) . . 

I jw I t9(t) dtds = -^'^ I w? ds -2m (57) 

Substituting this in the right side of (|56"D gives 



\ 2 f D yg(y) 2 f y ^f°tg(t) 2 dtdsdy = 

^g(D) 2 J D yg(y) 2 jj -^ds dy + \ f Q D g(y) 2 dy - ^Bl 



(5J 



2 ) jo jo 9 (s) 2 

Using this in (|56|) implies that the following inequality is sufficient for our result: 

\ 2 D 2 -l , , 2 f D , N2 r> 1 , , DgiDf , s 
g{D) 2 / yg{y) 2 / —ds dy - -^-J- > (59) 

1 Jo Jo 9\ s ) 4 

(for large D). 

We need only one additional fact about g(y); namely, that 1 > g{y) > e 2 . We 
already know that g' < giving the first inequality, and for the second, differentiate 

^— ' twice and use (^T|) to see that this expression, which is 1 at y = 0, then 



increases on [0, D\. 

Using this information it is seen that the double integral on the left of (??) does 
not tend to zero with D, and this proves (EO) for large D. This implies (|52"D and 



completes the proof of Lemma 16 



We also need a global bound on the solutions, for fixed k, d. 

Lemma 17 For given k and d, there is an M such that if (0, a, h) is a solution of 
(ED - (D mth > on [-d, d], then |0| + \<f/\ + \a\ + |a'| < M on [-d, d}. 

Proof: Letting ib = 2 — we see that ib'(±d) = and ib" > —k. This 

implies that ip' is bounded independent of which solution is being considered. Further, 
for any solution there is an x G (—d,d) with a(x ) = 0. If a'(x ) is bounded then 
we are done, so we can assume that there are solutions with o!(xq) arbitrarily large. 
Since a' is a minimum at Xq, this implies that for any A, and e, there is a A such that 
if a'(xo) > A, then the length of the interval in which \a\ < A is less than e. 
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Since ip f (±d) = we must have 



d 

(if) - tp 3 ) dx = ipa 2 dx. (60) 



d J-d 



Since max[_ djrf ] if>(x) = 1, and if)' > —k, there most be an e > such that for any 
solution of (]!]) — @ , if) > ~ in some interval fl of length e. For solutions with a'(xo) 
sufficiently large we must have \a\ large in at least half of Q, which means that the 



right side of ( 60]) can be arbitrarily large, while the left side is bounded by 2d. This 



contradiction proves Lemma 17 



We remark that M. K. Kwong gave a proof using Sturmian methods 
Completion of proof of Theorem || 



With M as in Lemma [17|, truncate (JJ) — (^]) by letting 

g(x) = min(x, M 2 ) 

and considering 

<f>" = K 2 (g(a 2 + 2 ) - 1)0 (61) 

and 

a" = g(cj>) 2 a. 

with boundary conditions © . With the notation of (||) let F = (F, G, H) = F(a, ft), 
where (3 = (/?, 7, 6). We use a truncation to ensure that whatever the initial condition 
at x = the solutions exist up to x = d. 

Since J > when a is close to 1, (when F(a, (0, 0, 80(a))) = 0), this solution 
t(a) = (0, 0, 5 (o;)) is non-degenerate and isolated and has Brouwer degree sgn J = 1 
(for the map F with fixed a). Similarly, if a is small and positive, then J < and 
t(a) has degree -1. Hence there is a change of degree along the branch of symmetric 
solutions between a small and a = 1. Hence, by a slight variant of Theorem 1.16 of 
Rabinowitz in [21], there is a connected set of solutions of F(a,/3) = branching off 



the symmetric solution (a,t(a)) at a point a G (0, 1) where J = 0, and this branch 
will either be unbounded in R A , or return to (0,0,0,0), or meet the symmetric branch 
again. 

Note that there are no solutions with h = 0, since then a would have to be of 
constant sign and we could not have f, <\> 2 a dx = 0. Also, solutions cannot bifurcate 
from the symmetric branch near a = or a = 1, because J 7^ there. 

The solutions of interest are positive. (That is, <f> is positive.) Solutions on the 
bifurcating branch of asymmetric solutions start off positive as the branch leaves the 
symmetric branch, from continuity. From (|l]) it follows that solutions can only fail 
to be positive along this branch if (f) — > (uniformly on [— d, d] ). Also, Lemma 
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shows the solutions must remain bounded, and in fact, remain within the truncated 
region where ([]]) — @ apply. If — > 0, then a — > h(x + c). 

It is conceivable that there are several points a where J = 0. However, the change 
of degree ensures that one branch does not return to the branch of symmetric solu- 
tions, and so, by Rabinowitz's global result, must tend to (0, h(x + c)) for some c / 0, 
that is, to a normal solution. In fact, by the symmetry of the problem, there will be 
two branches bifurcating from the same point, one tending to the unique bifurcation 
point from the normal solution with c > and the other to the symmetric reflec- 
tion of this solution around 0. (The uniqueness of this bifurcation point follows from 
Theorem |6|.) This completes the proof of Theorem |8| 

1.3 Conclusion: 

The initial motivation for this paper was Seydel's bifurcation diagram, Figure 1. 
Our goal was to prove that in some parameter range the problem could have as many 
as seven solutions (five essentially distinct). Unfortunately we have not achieved 
this goal. There are at least two features of Seydel's curve that seem important in 
obtaining such a proof. We would like to determine where on the symmetric branch 
the bifurcation to asymmetric solutions does exist, and we want to know the direction 
of bifurcation at this point. These remain challenges for future work. However, we 
have verified that for large nd the desired bifurcation from the symmetric branch 
occurs, and furthermore, there is a curve of asymmetric solutions going from the 
symmetric branch to the normal state. We have also shown how earlier results 
on bifurcation from the normal state can be obtained without the use of detailed 
asymptotics for the linear problem. 
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